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ABSTRACT 

In order to determine what processes govern the assembly history of galaxies with rotating disks, 
we examine the stellar mass Tully-Fisher relation over a wide range in redshift partitioned according 
to whether or not galaxies contain a prominent bulge. Using our earlier Keck spectroscopic sample, 
for which bulge/total parameters are available from analyses of HST images, we find that bulgeless 
disk galaxies with z > 0.8 present a significant offset from the local Tully-Fisher relation whereas, at 
all redshifts probed, those with significant bulges fall along the local relation. Our results support the 
suggestion that bulge growth may somehow expedite the maturing of disk galaxies onto the Tully- 
Fisher relation. We discuss a variety of physical hypotheses that may explain this result in the context 
of kinematic observations of star- forming galaxies at redshifts z = and z > 2. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: kinematics 
and dynamics — galaxies: spiral 



1. INTRODUCTION 

A major goal in galaxy evolution studies is to funda- 
mentally understand the ev olving dynam ical and mor- 
phological forms of galaxies (Roberts 1969). The favored 
method of tracking the assembly of stellar mass as a frac- 
tion of the total mass in rotationally-supported galax- 
ies is the redshift-dependent Tully-Fi sher (TF) re l ation, 
which was first explored at z ~ 1 by |Vogt et al.| ( |1996[ 
19971. Subsequent studies of the stellar mass (Af*)-TF 
relation at intermediate-to-high redshifts revealed scat- 
ters ~ 3x larger than that of the local relation (e.g ., 



Conselice et al.|[2005| iKassin et al.||2007] IVergani et al. 
2012 r This increased scatter was initially thought to 



represent a weaker coupling between stellar and dynam- 
ical mass, precluding detailed studies of the evolution in 
either slope or nor maliz ation. However, we showed in 
Miller et ail ([201 1[ ) and |Miller et ah] ( |2012| ) (hereafter, 
Mil and M12, respectively) with data of improved sig- 
nal/noise and refined modeling techniques, that the M„- 
TF relation is actually well-established at z ~ 1 with 
a scatter comparable to that seen in the local relation. 
Moreover, in M12, we demonstrated that the relation 
holds for most disk galaxies since z ~ 1.7, thereby pos- 
ing a challenge of how to explain the rapid evolution in 
kinematic behavior since z ~ 2 where star-forming galax- 
ies are morphologically-irregular and dispersion domi- 
nated (Forster Schreiber et al.||2006 Law et aL] |2007| 
Forster Schreiber et al.||2009p . ToTE e extent that aTTF 
relation can be examined at z ~ 2 ( Cresci et al.||2009"} 
Gnerucci et al.|20lT ) , a normalization increase of 0.4 dex 
is seen over ~1 Gyr to z ~1.5, in contrast to only 0.02 
±0.02 dex over the subsequent 9 Gyr. 

Since in M12 the TF scatter is observed to decline 
by 60% from z ~1.7 to z ~1, in this present paper we 
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seek to examine whether this might arise from physi- 
cal properties governing the evolution onto the TF re- 
lation. We target our attention on the morphological 
appearance of each galaxy, specifically the bulge-to-total 
ratio. Bulgeless disk s representing at least 15% of local 
galaxy populations (Kormendy et al. 2010) provide an 
interesting challenge for hierarchical ACDM galaxy for- 
mation (which leads to inevitable bulge-growth without 
substantial feedback: |Robertson et al.||2006| |Governatoj 
et al.||2010[ ). We test whether the high redshift M„-TF 
relation can be better understood when tracking the ma- 
ture, bulge-dominated population of galaxies separately 
from the evolving population of bulgeless systems expe- 
riencing a more secular formation p rocess. 



Throughout the paper we adopt a Chabrier ( 2003 ) ini- 
tial mass function and a J7a — 0.7, \ l m = U.3, H = 70 
km s~ x Mpc~ x cosmology. All magnitudes refer to those 
in the AB system (Oke 1974). 



2. DYNAMICAL DATA AND STELLAR MASSES 

The key measurements required to follow the evolving 
M*-TF relations are disk kinematics as parameterized 
through rotation curve model fits, and stellar mass esti- 
mates derived from multi-band photometric data. Our 
earlier papers (Mil, M12) describe the relevant data and 
their reduction in considerable detail so we provide only 
a brief summary here. 

Our spectroscopic sample was selected from Hubble 
Space Telescope (HST) Advanced Camera for Surveys 
(ACS) imaging data in various survey fields complete to 
an apparent magnitude of z=22.5 and is morphology in- 
clusive, containing irregular and merging systems as well 
as regular spirals with and without bulges. Keck spec- 
troscopy was undertaken for 236 galaxies with 0.2 < z < 
1.3 at a median spectral resolution of 30 km s _1 using 
the DEep Imaging Multi-Object Spectrograph(DEIMOS 
Faber et al. 2003) and, subsequently, 70 1.0 < z < 
1.7 galaxies were targeted at a median resolution of 
57 km s" 1 w ith the Low Reso lution Imaging Spectro 



graph (LRIS Oke et al. 1995) equipped with a red- 
sensitive CCD. An unique aspect of both spectroscopic 
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campaigns was the use of long exposure times (4-8 hrs) 
essential for tracking the rotation curves to the flattening 
radius (see Ml for details). Rotation curves were derived 
using various emission lines (Ha, [On], and [Oni] de- 
pending on the galaxy redshift. As discussed in Mil, wc 
account for position-dependent dispersion and emission 
brightness profile, convolved with the seeing, and adopt 
an arctangent function. We use inclination-corrected 
fiducial velocity measurements at 3.2 times the disk scale 
radius. The final sample for consideration here comprises 
171 galaxies for which rotation curves could be deter- 
mined (this is all galaxies except for spectrally compact 
or passive galaxies: see Mil & M12 for details). 

Stellar mass estimates are determined using a combi- 
nation of ground-based K-band infrared imaging, multi- 
band optical photometry, and spectroscopic redshift in- 
formation using the spectral energy distribution (SED) 
fitting technique first utilized by |Brinchmann fc Ellisj 
p000| . Measured magnitudes in multiple bands were 
applied u sing a Bayesian code ba sed on the precepts dis 



cussed in Kauffmann et al. ( 2003J), and later Bundy et al. 
(2005 1 . Using probability distribution functions that in- 



corporate uncertainties in the photometry, the stellar 
mass uncertainty is better than 0.2 dex for 83% of our 
sample. 

3. MORPHOLOGICAL DATA 

Our primary goal is to investigate the possible role 
bulge formation may play in the apparent rapid evolu- 
tion of the M*-TF relation from z ~ 2 to z— 1. W e facili- 
tate this investigation with GALFIT (Peng 2010). As we 
required disk scale lengths for earlier applications, the 
bulge-to-disk decomposition p rocedure desc ribed is sim- 
ilar to that in Mil, Ml 2, and|Miller| (|2012|) and so only 



briefly discuss the procedure here. 

We run GALFIT on each galaxy 1000 times, varying 
the initial parame ters in gaussia n distributions based on 
their SExtractor (Bertin 1996) values. For each object 
we attempt to fit a de Vaucouleurs bulge profile and an ex- 
ponential disk component, where the fit parameters are 
the center position, total magnitude m to t7 effective ra- 
dius R e (scale radius, r s , for an exponential disk), Sersic 
index n (fixed to n = 4 for deVaucouleurs and n = 1 for 
disk), axis ratio q, and position angle (j>. Where phys- 
ical bulge solutions are not found, we re-fit the galaxy 
with an index-varying single Sersic component (indices 
typically lie between 1 < n < 4). Such cases generally 
represent disk galaxies which are bulgeless and/or irreg- 
ular. Disk sizes, inclinations and position angles were 
taken from best-fit disk components if more than one 
component was fit. Final parameter uncertainties from 
the Monte Carlo distributions are better than 5% on av- 
erage, and we add these uncertainties in quadrature to 
the photometric errors from GALFIT. The scale radii, 
position angles and inclinations are typically measured 
better than 10%. Uncertainties are propagated through 
to TF parameters, resulting in larger errors for those 
galaxies which are difficult to constrain. 

In the DEIMOS sample, ^40% were adequately fit us- 
ing a 2-component decomposition, and ^60% benefitted 
from a single n- varying Sersic profile fit. In the LRIS 
sample, the relevant percentages were ~63% and ~37%, 
respectively. This serves as a good indication of the mor- 
phological distribution of our sample; less than half are 



well- formed spirals with a clear bulge (<|2j). Where HST 
data is available in multiple bands we compared gal- 
fit runs between bands to test for differences in the 
scale radius determination as a function of redshift. The 
scale radii are consistent among the bands indicating no 
significant redshift-dependent bias (less than 5% in the 
DEIMOS sample and <10% for that of LRIS). In order 
to maximize signal/noise we use the GALFIT results from 
the reddest available filter (F814W or F850LP). 

A crucial issue affecting classification at high redshift is 
the 'morphological k-correction' - the change in apparent 
morphology with increasing redshift following the drift 
blue ward in rest-frame wavelength. This is potentially 
troublesome for z > 1 where the the F814W and F850LP 
images sample the younger star-forming regions rather 
than the older, redder populations that dominate the 
stellar mass at lower redshift. As such, there is a danger 
of underestimating the bulge contribution. 

The HST near infrared Wide Field Camera 3 
(WFC3/IR) provides a F160W filter, which at 1 < z < 2 
provides rest-frame optical light and is therefore ideal 
for the bulge-to-disk decompositions we seek. While 
deep WFC3/IR F160W imagi ng from the CANDELS 
survey (Koekemoer et al. 20111 is available for one-fifth 



of our sample, the majority of the combined LRIS and 
DEIMOS samples are unfortunately in GOODS North 
(the WFC3/IR coverage of which will not complete for 
at least another year). However, for the purposes of this 
paper we seek only to demonstrate that use of the ACS 
data to classify the sample broadly into bulge-less and 
bulge-dominated subsets does not induce significant bi- 
ases. As we discuss below, we will split our overall sam- 
ple according to a dividing bulge-to-total ratio (B/T) 
= 0.1. With this division, we find, for the data with 
present WFC3/IR coverage, that morphological classifi- 
cations into these two categories are consistently made 
between the WFC3/IR and ACS data for 85% of the 
total sample. 

4. RESULTS 

We now examine the stellar mass Tully-Fisher (M*- 
TF) relation partitioned by morphology, in terms of the 
bulge/total ratio, B/T. To facilitate this investigation, 
we separate our sample according to the HST-derived 
GALFIT results into galaxies with prominent bulges and 
those without (bulgeless disks and irregulars) as de- 
scribed above. We plot the M*-TF relation in four red- 
shift bins (0.2 < z < 0.5, 0.5 < z < 0.8, 0.8 < z < 
1.2, 1.2 < z < 1.7) ensuring nearly equal sub-samples 
and look back time intervals (see Fig. [T|. Using the 
method described in Mil, we fit inverselinear regres- 
sions to each subsample and z-bin using a fixed slope (of 
3.70), the value of which was derived by fitting a free 
slope to the entire sample. In the two highest redshift 
bins (0.8 < z < 1.2, 1.2 < z < 1.7), we see bulgeless disks 
are significantly offset in the stellar mass (y-axis) normal- 
ization from that of the local relation by — 0.23±0.06 dex 
and —0.34 ± 0.07 dex, respectively. In contrast, disks 
with significant bulges do not deviate significantly from 
the local relation, nor in fact do bulgeless disks in the 
two lower redshift bins (Fig. §. The presence of a bulge 
appears to secure a disk galaxy on the M*-TF relation 
to within a scatter of 0.2 dex. 

The question then arises as to whether increased scat- 
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Fig. 1. — The M*-TF relations in four redshift bins. Galaxies with bulges (B/T > 0.1) are marked with red circles, and bulgeless/irregular 
galaxies (B/T < 0.1) are marked with navy boxes. The relation of bulge galaxies (red dashed lines, 1-cr) lie along the local relation in each 
2 bin, whereas the bulgleless relation (navy dotted lines, 1-cr) are offset from the local relation in the two highest 2-bins. To aid the eye, 
we plot the total sample with light grey error bars in each panel and include a solid grey line to mark the local relation. To understand 
the significance of the offsets, a Monte-Carlo 'bootstrap' analysis is conducted for each 2-bin, plotted in the inset panels along with the 
location and significance of the true offset. 



ter around the total Af*-TF relation from z ~1 to z ~ 1.7 
can be accounted for largely via the inclusion of bulgeless 
galaxies. In M12, the scatter from z ~ 1 to z ~ 1.7 in- 
creased up to 60%, whereas scatter across the 3 lower bins 
does not significantly evolve. Since the bulge-separated 
relations in the highest redshift bins have tighter rela- 
tions separately than when both samples are combined, 
it seems likely that increased scatter can be attributed 
to the zero-point shift of the bulgeless relation. The 
paucity of lower mass galaxies in the 0.8 < z < 1.2 
z-bin arising from the K-band magnitude limit applied 
for the DEIMOS sample likely complicates this inference 
(TF s catter increases to lower masses, e.g., |Begum et alT 
2008 1 . We note that the LRIS sample forming the basis 
ot the highest z-bin was not -FT-band limited. 

To allow for various differences between sample sizes 
and distributions, we quantify the offset significance for 
the bulgeless sub-sample using a Monte Carlo 'boot-strap' 



analysis. We fit two TF relations to randomly-selected 
subsamples (x 10,000, with replacement) of our full data 
sample (across all redshifts), ensuring subsamples of the 
same size (N) for each z-bin (N=21/22, 34/15, 26/20, 
21/12, for bulgeless/bulge-dominated sets respectively). 
With each pair of randomly selected subsamples, we mea- 
sure their relative normalization offset, and fit Gaus- 
sians to the resulting histograms of the offset distribu- 
tions (where the FWHM of each distribution is 0.082, 
0.099, 0.086, 0.112 dex for each bin in increasing red- 
shift). We also conduct an additional boot-strap analy- 
sis of N=10,000 where we select subsamples at random 
within redshift bins (and also with replacement). This 
latter method, which accounts for redshift-dcpcndence in 
the errors, results in broader normalization offset distri- 
butions (where FWHMs are 0.098, 0.110, 0.131, 0.132, 
respectively), and these distributions are plotted in the 
top left corner of each redshift-bin panel of Fig. [I] The 
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Fig. 2. — Evolution of the relative normalization of the M*-TF re- 
lation, defined as the difference from the local relation (zero-points 
from fixed-slope fits). Bulgeless galaxies (B/T<0.1) are denoted 
separately from the bulge sample (B/T>0.1). We indicate total 
relation scatter with the gold gradient fill; the offset normalization 
of the two subsamples occurs within the total scatter of the M»- 
TF relation. The bulgeless subsample becomes offset at high-z, 
but galaxies with bulges do not significantly deviate from the local 
relation (denoted by the solid black line, AAf*-offset = 0.0). 

true normalization offset observed in each bin relative to 
the bootstrapped distribution of the latter method yields 
offset significances of 1 .19c, 0.88cr, 3.14c, 3.07tr for each 
bin in increasing redshift. Together this translates to a 
confidence interval of greater than 99.8% that the rela- 
tion of the bulgeless galaxies are offset at high-z due to a 
genuine effect rather than random error or scatter. The 
slight decline in significance in the highest bin is due to 
the reduced number of galaxies with bulges in that bin, 
even though in real terms the offset of the bulgeless re- 
lation is greatest in the highest z-bin (— 0.34±0.07). 

5. DISCUSSION 

In light of our results, we explore explanations of the 
offset co-evolution of stellar mass and total mass assem- 
bly in bulgeless galaxies from those which experienced 
earlier bulge growth. We first consider in § 5.1 a favored 
picture in the literature of z ~ 2 studies, and then turn 
in §5.2 and §5.3 to what may be a more self-consistent 
and simple picture for our results. 



5.1. Clump formation and migration 
2, low ratios of rotati on-to-dispersion velocity 



At z 

support in disks {V/cr < 2 — 5, |Forster Schreiber et al 
2006 2009) have been interpreted in a simple picture 
whereby disks fragment and the resulting clumps migrate 
to the center s of ga laxies to for m bulg es ( |Noguchi||1999| 
Immeli et al ||2004| [Romeo et al.|2010 |. 

Importantly, the velocity dispersion measured from 
such studies comes from the ionized gas. The disper- 
sion in the emission lines could be dominated by energy 
or momentum driven stellar feedback, rather than dis- 
persion from dynamical pressure (as would be traced by 
the stars and cold gas). In our sample, while there is a 
spread in V/cr values from 2-15 across all redshift bins 
(90% never rise above V/a=10), no clear trend in a it- 
self exists with respect to morphology or redshift. These 
values are more similar to the locally observed spread 
in V/cr of ionized gas than those found in massive, star- 



Furthermore, the clump migration picture does not ex- 
plain why bulgeless galaxies arrive on the M^-TF relation 
by intermediate redshift without forming bulges. In the 
clump migration picture, gas-rich bulgeless and barless 
disks would continue to suffer gravitational instabilities 
and clump formation until a bulge had formed and sup- 
posedly stabilized the disk. Some other mechanism is 
needed to explain the stabilization of bulgeless disks that 
have yet to form a substantial bulge/bar. 

Additionally, central stellar velocity dispersions can 
support an increased stellar mass in the form of an accu- 
mulating bulge without significantly changing the rota- 
tional velocity of the surrounding disk. This would shift 
galaxies above the TF relation as their bulges grow with- 
out increasing rotational support, which is not supported 
by our results. 

5.2. An underestimation of gas masses in high-z 
bulgeless disks 

A simpler explanation for our results is that bulgeless 
galaxies have higher gas fractions in their disks as com- 
pared to the rest of the sample. By adding total gas 
masses to our stellar masses, we may find the baryonic 
Tully-Fisher is universal at all redshifts in rotationally- 
supported galaxies. To explore this possibility without 
direct gas mass measurements, we conduct the following 
exercise. 

We compare the gas estimates from the empirically- 
based analytical method from Mil (Method 1), to the 
estimates based on the Kennicutt-Schmidt (K-S) rela- 
tion (Method 2). Described in detail in Mil, Method 1 
uses local gas fractions based on a galaxy's stellar mass 
and self-consistently adds the integrated specific star- 
formation rate from the galaxy's redshift (accounting for 
a gas recycling fraction based on the assumed Chabrier 
IMF). In Method 2, the K-S gas masses are determined 
from the rest-frame £?-magnitude surface brightnesses, 
which are used to estimate a SFR density and thus a gas 
mass according to the inverted K-S relation ( |Kennicutt 
1998]) . 

This comparison reveals that the high redshift bulge- 
less disks have a significantly diminished B-band magni- 
tude per area compared to the rest of the sample. That 
is to say Method 1 results in 1.11 times more gas than 
Method 2 with a 10% scatter in bulgeless disks, whereas 
this factor is 1.43 in the rest of the sample. If we in- 
crease the bulgeless disk gas mass estimates 30% so that 
they align with the rest of the sample, then a universal 
baryonic TF relation is restored. 

Physically, a correction of this nature to our gas mass 
estimates suggests less metal-enriched gas in bulgeless 
disks, which is less efficient at forming stars. The rel- 
ative youth of the bulgeless disks may be due to the 
lack of enriching outflows re-condensing at z ~ 1 (but 
do by z ~ 0). A lengthened 'fountain' duty-cycle could 
reflect a shallower gravitational p otential and lower star- 
formation surface densities ( i.e., Oppenheimer & Dave 



2008 Finlator & Dave 2008 etc.), where outflows are 



slower to re-condense and metals more likely to escape 
in supernova-driven winds. 

Also the simultaneous growth of supermassive black 
holes (SMBHs) at the centers of galaxies with their bulge 
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mass is well known (Ferrarese & Merritt 20001. SMBHs 
are likely fed by gas that has sunk to the ;enters of galax- 
ies via disk instabilities and/or mergers, a process which 
would similarly grow a bulge. It is then unsurprising that 
galaxies with no central bulge may have a much larger 
fraction of their gas still in their disks. 

While the predictions of this exercise await testing via 
direct gas mass observations, it does suggest interesting 
implications for a universal baryonic-TF rela tion (so far 
confirmed only locally, e.g. McGaugh 2012). Looking 
towards the future, these predictions can be tested by 
determining the molecular and neutral hydrogen com- 
ponents of these galaxies, via facilities such as the At- 
acama Large Millimeter/submillimeter Array (ALMA), 
or future radio facilities (e.g., MeerKAT, or ultimately 
the Square Kilometer Array, SKA). 



(appropriate for our study since the role of environment 
would be subtle in our field sample). These models are 
ultimately driven by the cosmic decline in accretion 
rate, shutting down assembly of massive galaxies first 
by quickly consuming their reservoirs. A combination of 
the evolving ultraviolet background with photo-ionizing 
radiation from the first stars could create a transition 
mass, above which the cooling efficiency is relatively 
higher, and below which a lack of self-shielding keeps 
smaller, thinner disks from remaining neutral (keeping 
molecular gas collapsing to form GMCs). Galaxies 
with bulges maintaining thicker disks in thicker, steeper 
potential-wells could self-shield, and thus form stars 
more efficiently than thin, bulgeless disks being ad- 
versely affected by photo-ionizing radiation in shallower 
potential- wells. 



5.3. A z-dependent transition mass 

Locally, low-mass bulgeless galaxies (M* < 1O 9 M ) 
tend to fall bel ow the extrapola ted M*-TF relatio n from 
M* > 10 9 M Q {Matthews et al.||l998| |Stark et al.||2009|). 



Since this offset is similar to that seen in our high-redshirt 
bulgeless galaxies, it suggests that probing further down 
the stellar-mass function at 0.2 < z < 0.8 may uncover 
a similar transition to an offset of bulgeless disks at an 
intermediate mass between that observed at z ~ and 
atz > 1. 

The physical significance of an evolving transition 
mass for bulgeless galaxies could be understo od via 
the 'downsizing' concept (e.g. Cowie et al. 1996) or an 
evolving 'mass floor ' in galaxy formation theory (e.g. 



Bouche et al. 2010). These models seek to explain 
why more massive galaxies formed earlier and faster 
than lower mass galaxies, regardless of environment 



In attempting to better understand drivers of disk as- 
sembly from our results, we note a number of tensions 
regarding the picture where disks settle from the migra- 
tion of large clumps into central bulges. Rather, a more 
self-consistent picture could be provided by an univer- 
sal baryonic-TF relation, where better accounting of gas 
in high-z disks could explain the offsets seen in tracking 
stellar mass with the total rotational support in galaxies. 
This may also predict a redshift-dependent transition 
mass which lowers with the age of the universe, below 
which bulgeless disks assemble their mass offset from the 
locally-defined M*-TF relation, but not the baryonic-TF 
relation. 
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